Atomistic spin-model based on a new spin-cluster expansion technique: Application to 

the IrMn 3 /Co interface 
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In order to derive tensorial exchange interactions and local magnetic anisotropies in itinerant 
magnetic systems, an approach combining the Spin-Cluster Expansion with the Relativistic Disor- 
dered Local Moment scheme is introduced. The theoretical background and computational aspects 
of the method are described in detail. The exchange interactions and site resolved anisotropy con- 
tributions for the IrMn3/Co(lll) interface, a prototype for an exchange bias system, are calculated 
including a large number of magnetic sites from both the antiferromagnet and ferromagnet. Our 
calculations reveal that the coupling between the two subsystems is fairly limited to the vicinity of 
the interface. The magnetic anisotropy of the interface system is discussed, including effects of the 
Dzyaloshinskii-Moriya interactions that appear due to symmetry breaking at the interface. 
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I. INTRODUCTION 

Spin models form an important tool in the understand- 
ing of magnetic systems, particularly of multilayers sys- 
tems such as bilayers consisting of ferromagnetic and an- 
tifcrromagnctic materials, because such models allow the 
investigation of magnetic configurations and dynamics at 
the atomic scale. Bilayers of ferromagnetic and antifcr- 
romagnetic materials are important parts of magnetore- 
sistive sensors, used, e.g., in hard disc read heads. The 
role of the antiferromagnetic component is to pin the fer- 
romagnetic layer by the so-called exchange bias effect)^ 
leading to a shift of the hysteresis loop of the ferromag- 
net. A number of calculations based on spin models have 
revealed that the anisotropy of the antiferromagnet and 
the strength of the interface coupling between antiferro- 
magnet and ferromagnet are crucial for the strength and 
stability of exchange bias 3 - - — and therefore for the func- 
tioning of the device. 

IrMn is the most common antiferromagnet industrially 
because of the large exchange bias and thermal stability 
that can be obtained with this material. Despite the 
successful application of this material in exchange bias 
bilayers, there is still no complete microscopic under- 
standing of the most important parameters, the antifer- 
romagnetic anisotropy and the interface coupling. These 
properties have recently been inferred by a careful ex- 
amination of the temperature behaviour of the exchange 
bias of polygranular thin film bilayers^. In order to evalu- 
ate the antiferromagnet anisotropy, the magnitude of the 
exchange bias is fitted to a simple model that depends 
only upon the grain size distribution, yielding effective 
anisotropy and coupling parameters. The origin of the 
effective anisotropy and interface coupling, due to the 
microscopic physical properties of the interface, remains 
however inaccessible to direct experimental study. 

On the theoretical side, only few first principles calcu- 



lations on relativistic effects in Mn-based antifcrromag- 
nets have been performed. Regarding the anisotropy, 
Umetsu et al& calculated the magnetic anisotropy energy 
(MAE) of Llo-type MnTM (TM=Ni, Pd, Pt, Rh and Ir) 
compounds and revealed that Mnlr had the highest MAE 
with a value of -7.05 meV/unit cell. Furthermore, a gi- 
ant second order magnetic anisotropy for LI 2 IrMn3 was 
reported, leading to energy barriers of the order of K e & 
= 10.4 meV/unit cell for rotation of the Tl ground state 
spin-structure around the (111) axis£ To the best of our 
knowledge no first principles calculations of the interface 
exchange coupling between an antiferromagnet and a fer- 
romagnet have been published to date. 

The organization of the paper is as follows: in the first 
part we introduce the concept of the Spin-Cluster Ex- 
pansion (SCE) merged with the Relativistic Disordered 
Local Moment (RDLM) scheme and discuss the details 
of the first-principles calculations. Then we present our 
results for the magnetic moments, the exchange coupling 
and the magnetic anisotropy in LI2 IrMn3, fee Co bulk, 
and for the IrMn3/Co(lll) interface. Special attention 
is paid to the effect of Dzyaloshinskii-Moriya interactions 
at the interface. 



II. THEORY 

A. The Spin-Cluster Expansion 

Mapping of the energy of an itinerant magnetic system 
onto a spin-Hamiltonian presents a longstanding problem 
in computational solid state physics. The most widely 
used approach relies on an adiabatic decoupling of the 
longitudinal and transverse spin- fluctuations, determined 
by the fast motion of mobile electrons, r e ~ 10~ 15 s, 
and the slow motion of spins, T s f ~ 10 s. Under 
such circumstances, Time-Dependent Density Functional 
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Theory can be approximated by a Landau-Lifshitz equa- 
tion of motion for the classical spin variables, i.e., the 
orientation- field, e(f , t) , as combined with Spin-Density 
Functional Theory (SDFT) within the Local Spin Den- 
sity Approximation (LSDA)J^ Further simplification of 
the theory is represented by the Rigid Spin Approxima- 
tion, in which the orientation of the spins is allowed to 
change from atomic cell to atomic cell. A system of N 
spin-moments is thus characterized by the set of orien- 
tations (spin-configurations) {e} = {ei, . . . , ejv}, with 
|ej| = 1. The energy (grand-potential) of the system, 
f2({e}), determined within LSDA then defines a classical 
Hamiltonian that can be used in spin-dynamics simula- 
tions. 

In order to ensure that the process is tractable, a suit- 
able parametrization must be applied to fi({e}). The 
simplest reliable approximation, an (extended) Heisen- 
berg model, is of second order: 

JV N 

f2({e}) = fi + 2_j ^ K ^i + 2 E ' "' <•/''/ ' ( 1 ) 

(MO 

where Kj and Jjj are the second-order anisotropy matri- 
ces (which are traceless and symmetric) and the tenso- 
rial exchange interactions, respectively. The Relativis- 
tic Torque Method (RTM), which was introduced in 
Refs. ITTlfl2l . allows the evaluation of these parameters us- 
ing relativistic electronic structure calculations. By using 
RTM, one is required to choose a reference magnetic con- 
figuration; usually the ferromagnetic state is taken. For 
complex magnetic structures, however, the use of such a 
reference state might lead to incorrect results concerning, 
e.g., the magnetic ground state of the system Another 
drawback of RTM is that for systems with reduced point 
group symmetry the derivation of the anisotropy matrix, 
Kj, might be quite a cumbersome task. More impor- 
tantly, RTM is not directly applicable to the calculation 
of more complicated two-spin interactions, such as bi- 
quadratic coupling, or multispin interactions, known to 
be present in magnetic systems at the nano-scalei 13 ' 14 

The Spin-Cluster Expansion (SCE) developed by 
Drautz and Fahnl o 15 i 16 in principle solves the prob- 
lems mentioned above, in that it provides a system- 
atic parametrization of the adiabatic magnetic energy 
of general classical spin-systems. In the following we 
briefly outline the method. Considering clusters of n sites 
(0 < n < A), C„ C I N = {1, 2, . . . , A}, the functions 

"v. : ' ({*» = uJ N - n) n II ^ . ( 2 ) 

L ik = (£ ik ,m ik ) ^ (0,0) 

form an orthonormal and complete basis set of the func- 
tion space over {e}. Here Y^(e) denote, for convenience, 
real spherical harmonics. The grand-potential of the sys- 



tem can, therefore, be readily expanded as 

= E E J ' : ' ' II Y ^ & k ) , 

C n L il ,...,L i „/(0,0) i t ec„ 

(3) 

where the multispin interactions, J^ 11 ' "' are defined 
as 

' = / d 2 e h ...d 2 e in (ft) eV .. e%> I] Y ^ k 

i fc ec„ 

(4) 

In the above equation, the key quantities of the theory, 
namely, the restricted averages of the grand-potential 
are: 

We V ..e in = TT^T / d 2N - 2n e lNXCn Sl({e}) , (5) 

1 (47T) J 

i.e., by fixing the spin-orientation at the sites of cluster 
C n , an average must be taken on all spins outside of the 
cluster, i € In \ C n . Restricting ourselves to one-site and 
two-site interactions, Eq. ^ reduces to 

nm) ~ n + J2 E ( 6 ) 

i L#(0,0) 

+ It, e e wylwyv®), 

with 

n = (n) , (7) 

Ji=j d 2 e t (n),Y L (e, t ) (8) 

and 

4 L ' = J d 2 e t J d 2 e 3 Y L (e,-) Y L i % ) . (9) 

The spin-Hamiltonian, Eq. (Q]), can easily be related to 
Eq. ^ through the relationships, 

2 

ei K 1 e l = ^2 J> 2 ' m ' Y 2 , m (ei) + const. (10) 

m=-2 

and 

Ci Jij Cj > J^j ^l,m *l,m' ) • 

m,m , — —l 

(11) 

The evaluation of Eq. ([5]) obviously poses a heavy 
numerical task since it involves integration of the first- 
principles grand-potential over N — n spin-variables. In 
particular, for infinite systems like bulk magnets or mag- 
netic films, it seems extremely demanding to ensure the 
completion of this task. A suitable method is therefore 
needed that can evaluate the restricted averages in a fea- 
sible manner. 
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B. The Relativistic Disordered Local Moment 
scheme 

The Disordered Local Moment (DLM) scheme^ is 
an extension to conventional SDFT to include trans- 
verse spin-fluctuations. In the DLM scheme the mag- 
netic excitations are modelled by associating local spin- 
polarization axes with all lattice sites and the orienta- 
tions {e} vary very slowly on the time-scale of the elec- 
tronic motion. These "local moment" degrees of free- 
dom produce local magnetic fields centered at the lat- 
tice sites which in turn affect the electronic motions 
and are self-consistcntly maintained by them. The the- 
ory of Disordered Local Moments (DLM) in conjunc- 
tion with the Korringa-Kohn-Rostoker Coherent Poten- 
tial Approximation (KKR-CPA) was proposed twenty- 
five years ago by Gyorffy et alJ^ and has recently been 
generalized by Staunton et al.— ~— to include relativistic 
effects. Here we give a short summary of the method as 
applied to the consideration of the paramagnetic state. 

For spherically symmetric potentials and exchange 
fields, the local orientation of the spin-magnetization is 
accounted for by the similarity transformation of the 
single-site f-matrix, 



U {e i )=D{e i )t i {e z )D{e i y 



(12) 



where for a given energy (not labelled explicitly) i 4 (e 2 ) 
stands for the i-matrix with exchange field pointing along 
the local z axis and D_ (e*;) is a (block- wise) projective rep- 
resentation of the SO (3) transformation that rotates the 
z axis onto e%. Note that underlining indicates matrices 
in the (k, /j,) angular-momentum representation. In the 
spirit of the single-site coherent-potential approximation 
(CPA) an orientation independent reference medium is 
introduced in terms of the effective i-matrices, c , and 
the corresponding scattering path operator (SPO) ma- 
trix, 



l c = ( f 



(13) 



In the above equation, double underlines denote matri- 
ces in site-angular momentum space; t is diagonal in 
site indices while G Q stands for the matrix of structure 
constants. The effective ^-matrices are determined by the 
condition that the excess scattering matrices, 



(14) 



vanish when averaged over all possible orientations with 
an equal probability in the paramagnetic state, 
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2Q (Si) d 2 e t = 



(15) 



The above condition can also be reformulated in terms 
of the site-diagonal SPO matrices, 



(16) 



where the so-called restricted average of the site-diagonal 
SPO matrix can be expressed as, 



with 



Di(?i) = L + Ki (3s) i», 



t,c } — %i,c\ 



(17) 



(18) 



As discussed in the previous section, we limit our con- 
sideration by neglecting longitudinal spin-fluctuations; 
i.e. we suppose that the LSDA potentials and exchange 
fields are independent of the orientational state {<?}. The 
magnetic force theorem^ then implies that the LSDA to- 
tal energy can be replaced by the single-particle energies 
(band energy) and the grand-potential at zero tempera- 
ture can be expressed as 



(M) = E LSDA ({4) - e F N ({e}) 



deN(e;{e\) . 



(19) 



Here Ef is the Fermi energy and N(£;{e\) denotes the 
integrated density of states. Using Lloyd's formula^ and 
taking the coherent medium as reference, Eq. (fl~9| can be 
recast as 

n ({el) = n c - i Im y de IndetZ). (3) 

i 

fe=l ii^i2^...^ik 

Tr {X n (S il )r cMi2 X i2 (e l2 ) ...X ik (3 J r^.J , (20) 
with 

Q c = -J F d£N {£)-^-Im J F ds lndctx c (e) . (21) 

where No (e) represents the integrated density of states 
of the free electrons. 

Eq. (|20[) is, in principle, valid for reference systems 
fixed by any set of t- matrices, t .. The particular choice 
of the RDLM effective medium, defined by the condi- 
tion Eq. (|15p , is very useful in finding a suitable approx- 
imation for the restricted averages in Eq. ([5]) and conse- 
quently for the spin-cluster interactions in Eq. (jj). More 
specifically, when taking averages we neglect the so-called 
back-scattering contributions in the third term of the 
right-hand side of Eq. (|20[) . i.e., those for which any site 
outside of the chosen cluster, i; e Ijv \ C n (I = 1, . . . , k), 
occurs at least twice. Under this restriction, which is 
consistent with the single-site CPA^ we obtain for the 
constant term in Eq. ([7|), 



(22) 
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The one-site restricted average of the grand potential can 
be expressed as, 

1 f F 

(fi)- =n c Im / de In det (e*) 

- \ E Im / d£ ^J Se i lndct ^- fe) > 



(23) 



from which one obtains for the one-site expansion coeffi- 
cient, 

= -~Im / J d 2 e t lndet£>i (e 4 ) F £ (ei) . (24) 
From the two-site (i =/= j) restricted averages, 

= 

f f ep 

Im / de (In det D t (e*) + In det D_ (ej)) 

Im / de ^2 — / ^ 2e ' lndeti2;(e;) 

1 y 00 ^ 1 f EF r 



7T ' — ' k 

fe=i 



(25) 



and using Eq. ([S]) the following closed expression for the 
two-spin interactions can be deduced, 

J^ L ' =--Im J de J J d 2 ei d 2 e 3 

Trln (/ - X, (Si) T ctj X 3 %) r^) F L (e 4 ) Y L , (e}) . 

(26) 

In the Appendix we derive the non-relativistic limit of 
the above expression. 

It is also straightforward to show that a general multi- 
spin interaction as defined in Eq. ^ is given within the 
RDLM scheme as 



r L H „„M 



[■■ [ d 2 e il ...d 2 e in f[Y Lii (e H ) 
J J i=i 
Tr (X h (e h ) r clil2 X l2 (e h )... X ik (e lk ) r clfcn ) , (27) 

where each term on the right-hand side must contain all 
the sites in cluster C n . 



the IrMri3/Co(lll) interface in terms of the Screened 
Korringa-Kohn-Rostoker (SKKR) metho d 23 ' 24 . Because 
the magnetic structure of the IrMn3/Co(lll) interface 
was not known a priori, we employed in these calcula- 
tions the scalar-relativistic DLM approach that is repre- 
sentative of the paramagnetic state of the system. Here 
we solved the CPA condition, Eq. (|A12[) . with a relative 
accuracy of 10 -5 for the effective t-matrices. The local 
spin-density approximation as parametrized by Vosko ct 
ali2£ was applied, the effective potentials and fields were 
treated within the atomic sphere approximation with an 
angular momentum cut-off of £ ma x = 2. The energy in- 
tegrations were performed by sampling 12 points on a 
semi-circular path in the upper complex semi-plane and 
84 points were selected in the 2D Brillouin zone, i.e., in 
one fourth of the fcc(lll) surface Brillouin zone, for the 
necessary fc-integrations. In all cases we used a parent fee 
lattice structure with the lattice constant of LI2 LM113, 
a = 3.785 A 26 ' 27 , with no attempt at geometric relax- 
ation. 

To ensure that both the IrMn3 and Co components 
shared the same two-dimensional translational periodic- 
ity, which is necessary within the layered SKKR method 
for an interface, the calculations for Co were carried out 
with four atoms per cubic unit cell according to the fee 
lattice structure. The interface calculations were per- 
formed for a symmetric system consisting of two 6 mono- 
layer (ML) thick RM113 layers sandwiching a 12 ML Co 
layer from both sides. Considering four atoms in a unit 
cell for each layer (three Mn and one Ir atoms in each 
IrMn3 layer and four Co atoms in each Co layer), the 
total 2D unit cell of the interface system contained 96 
atoms. In addition, to ensure smooth boundary condi- 
tions, the entire system was sandwiched between two per- 
fect semi-infinite bulk IrMn3 systems as required by the 
SKKR method^s. 

In order to evaluate the integrals over the orientations 
in Eqs. and (f2"6"|) we used the efficient Lebedev- 

Laikov schemed The energy integration in Eq. (|2~4"|) was 
performed by introducing a damping through the Fermi 
function with T = 50K. The corresponding integral can 
then be transformed into a sum over the fcrmionic Mat- 
subara poles and the real part instead of the imaginary 
part of the logarithm of complex numbers has to be con- 
sidered, thereby avoiding the well-known phase problem 
associated with Lloyd's formula^ In Eq. (|2"6"|) the above 
mentioned contour integration can safely be used, since 
the expression, Tr In (7 — A), in the integrand can be ex- 
panded into a rapidly converging power scries. 



III. APPLICATION TO THE IRMN 3 /CO(lll) 
INTERFACE 

A. Computational details 

We performed self-consistent calculations for the or- 
dered bulk LI2 IrMn3 alloy, for fee bulk Co and for 



B. Magnetic moments in the DLM state 

The self-consistent field calculations resulted in a local 
magnetic moment of 2.20 /j,b for the Mn atoms in LI2 
IrMn 3 . This value is about 15 % less than the value we 
obtained for the the spin-moment in a rclativistic cal- 
culation for the ordered triangular (Tl) stated Such a 
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softening of the magnetic moments is expected in the 
paramagnetic DLM state due to the vanishing Weiss- 
field acting on the spins. This reduction of the local 
moment is not necessarily accompanied by a softening of 
the exchange coupling, however, as is demonstrated be- 
low. The calculated local magnetic moment of the fee 
bulk Co in the DLM state is 1.62^s, in agreement with 
the experimentally and theoretically obtained values for 
ferromagnetic bulk Co. 

In Fig.[T]the magnetic moments for the IrMn 3 /Co(lll) 
interface system are shown as calculated in the DLM 
state. Note that only one half of the symmetric sys- 
tem is shown. Clearly, only moments in close vicinity of 
the interface are substantially different from their bulk 
values: the Mn moments increase to 2.54^b and the Co 
moments reduce to 1.42/z.b- 
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A very good agreement between the SCE-RDLM and the 
RTM can also be seen in this case. Importantly, both 
systems are characterized by short-range interactions: 
the Jtj's practically vanish beyond the fourth neighbour 
shell. 
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FIG. 1: (Colour online) Calculated magnetic moments of Mn 
and Co atoms for the IrMn3/Co(lll) interface. The IrMn3 
and Co sides of the interface, marked by the vertical dotted 
line, correspond to layers labelled by 1 to 6 and 7 to 12, 
respectively. 



C. Isotropic exchange interactions 

Recently, we investigated the magnetic ground-states 
of the ordered IrMn 3 alloy^ In agreement with experi- 
ment^ and other theoretical work2£, we found that the 
ground-state spin-structure of the Mn sublattices is a 
frustrated 120° N eel- type state in the (111) planes of 
the fee lattice, called the Tl state. This frustration is 
the consequence of the strong antifcrromagnctic exchange 
coupling between the Mn atoms. In Fig. [5] we compare 
the isotropic exchange interactions, = (1/3) TrJ.^ , 
obtained from the newly developed SCE-RDLM scheme 
with those from the RTM^i Clearly the two sets of pa- 
rameters (circles and squares) are in very good quanti- 
tative agreement. This indicates that the softening of 
the magnetic moments in the DLM state do not imply a 
softening of the exchange couplings with respect to the 
ordered state. In contrast to IrMn 3 , the exchange inter- 
actions of bulk Co are primarily of ferromagnetic nature. 



FIG. 2: (Colour online) Calculated isotropic exchange inter- 
actions, Jij, for bulk LI2 IrMn3 (upper panel) and fee bulk Co 
(lower panel). Squares and circles refer to the Jij's obtained 
from the DLM states by using the SCE-RDLM method and 
by using the RTM,— using the ordered ground state (the Tl 
state for IrMn3 and the ferromagnetic state for Co) respec- 
tively. The interactions are displayed as a function of the 
distance, Rij. 

For the 6ML IrMn 3 /12ML C0/6ML IrMn 3 interface 
system we considered all the pairs of magnetic atoms sep- 
arated by a distance less than 10.68A. This implied the 
calculation of a total number of more than 13,000 ten- 
sorial interactions. Inspection of the calculated isotropic 
exchange interactions reveals that the Mn-Mn and Co- 
Co interactions are spread around their corresponding 
bulk values due to symmetry breaking at the interface: 
the nearest neighbor (NN) Mn-Mn interactions are some- 
what increased while most of the NN Co- Co interactions 
are decreased. This decrease is of up to 20% near the 
interface. This feature clearly correlates with the change 
of the local magnetic moments caused by the presence 
of the interface depicted in Fig. [TJ The Mn-Co interac- 
tions across the interface are quite small in magnitude (< 
10 meV) and mostly antifcrromagnctic: = — 6. 15 me V 
for nearest neighbours, while Jij = —4.86 meV and 
—8.40 meV for next nearest neighbours. 
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In order to visualize the isotropic exchange interactions 
for the interface system, we defined the following effective 
interaction parameters, 



T> 



'I -J 



(28) 



where p and q denote layers, while i and j stand for sites 
within the corresponding layers. For this calculation we 
choose a magnetic state, {e p i}, corresponding to the Tl 
structure for the IrMn3 side and to an FM state perpen- 
dicular to the interface within the Co layers. This con- 
figuration was chosen because of the antiferromagnetic 
nature of the interface coupling combined with the large 
Mn anisotropy, which constrains the Tl structure within 
the interface plane. Clearly, by translational and rota- 
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FIG. 3: Calculated effective exchange interactions, J p i, for 
the IrMns/Co interface as defined in Eq. (|28j|. The labelling 
of layers is identical to that in Fig. [1] 

tional symmetry of this spin-structure, J P i turned out to 
be the same for all the Mn sites in a given layer p of the 
IrMn3 region, while in each layer of the Co region, in 
principle, two different values for J p i may be obtained. 
In particular for the Co layer adjacent to IrMn3, these 
two different values are associated with Co sites having 
two and three NN Mn atoms. 

Fig. [3] clearly shows that the effective exchange inter- 
actions quickly approach their bulk value with distance 
from the interface. This follows from the fact that only 
interactions between atoms in close proximity to the in- 
terface deviate from their counterparts in the bulk; those 
more than 2 atomic planes away from the interface are 
practically bulk-like. Such behavior is to be expected be- 
cause of the very narrow range of exchange interactions 
in the bulk systems. Since for the chosen magnetic states 
the magnetic orientations of the Mn and Co atoms are 
perpendicular to each other, the otherwise weak Mn-Co 
interactions are excluded from the sum in Eq. (|28[) and 
the the effective interactions are significantly decreased 
near the interface. The different values of J P i for the 
non-equivalent Co sites near to the interface can also be 
seen in the figure. 



D. Magnetic anisotropies 

The most remarkable observation of our previous 
study^ was the surprisingly strong, second order MA of 
IrMn3 resulting from the fact that the cubic symmetry 
is locally broken for each of the three Mn sub-lattices. 
Rotating the Tl state around the axis normal to the 
(111) plane, the energy of the system can be described 
as E{ip) = Eq + iieff sin 2 ip. By calculating E(tp) in the 
Tl spin-state, a value of K c g=10A meV was deduced £ 
We repeated this calculation by using the self-consistent 
field potentials from the DLM state and obtained a value 
of K e s = 7.67 meV. Although this value is about 30% 
less than the one we obtained in our previous calculation 
using the Tl ground-state, it is still consistent with the 
observation of a very large magnetic anisotropy. 

The SCE-RDLM scheme provides a unique opportu- 
nity to calculate contributions of the MAE related to 
the one-site and the two-site exchange anisotropy terms, 
which is an important issue for spin-model simulations 
because it affects the scaling of the anisotropy with mag- 
netisation. First we note that, by symmetry, in the frame 
of reference with z axis along the (111) direction, the one- 
site anisotropy matrices of the three Mn sublattices in the 
LI2 phase can be expressed as, 
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(31) 



Our calculations based on Eqs. (fTUj) and (j2"4"]l confirmed 
with high accuracy the above non-trivial forms of the 
one-site anisotropy matrices, with a value of K = 0.82 
meV. Note that the total one-site contribution to K e s 
is In addition, from the tensorial exchange in- 

teractions, Eqs. ([2]) and ([TTj). we calculated a two-site 
anisotropy contribution of 6.73 meV to the MAE. Thus, 
the spin-Hamiltonian derived from the SCE-RDLM ap- 
proach yielded K c ^ = 8.37 meV, in good agreement with 
A' ff = 7.67 meV we calculated directly from the band- 
energy. The difference between the two values should be 
related to the fact that the direct calculation of the MAE 
is not restricted to second order interactions as is the case 
for the spin-Hamiltonian. Due to cubic symmetry, for fee 
bulk Co we calculated negligible MAE, \K e s\ < 1 fieV. 

Next we investigated the magnetic anisotropy of the 
spin-structure described in the context of Eq. (|28|) and 
Fig. El for the interface system. Of particular interest is 
the variation of the energy, Eq. (TT|), during a magneti- 
zation reversal process. For transparency, we considered 
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the reversal of the Co moments from z to —z by keep- 
ing the relative orientations of any two moments in the 
system fixed. Theoretically, this is equivalent to rotating 
all the spins simultaneously around an axis in the (111) 
plane by an angle of < (p < tt. Clearly from Eq. (fTJ the 
energy, E((p), can be recast into layer-wise contributions, 
E p (ip), such that 



mcV, KH 



=0.54 mcV and K P J S ~- 



(32) 



where N is the number of cells in a layer and E p (ip) com- 
prises the contributions of three Mn atoms and four Co 
atoms in the IrMn 3 and Co parts of the interface system. 
Moreover, E p (ip) can be decomposed into contributions 
arising from the one- and two-site terms of Eq. (jTJ) . 
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FIG. 4: (Colour online) Calculated effective anisotropy pa- 
rameters (triangles), K% s , see Eq. (|32[) , for the IrMns/Co in- 
terface. Also shown are the corresponding one-site (K' t 
circles) and two-site (K^ s , squares) contributions. 



p,os 

ff 



Note that in these calculations only the symmetric part 
of the tensorial exchange matrices, Jfj = | (jy + J^), 
were considered. Our numerical results confirmed that 
the following forms of E p (tp), which describe the bulk 
anisotropy energies^ give a good description for the lay- 
ered system as well, 



E p (fp) = E p (0) + ^(2 + 



sin 2 tp 



' 8 

— 2\/2sin<^(l — cos</?) 

for the IrMn3 layers and 

E p {<p) =E p {0)+K* s sm 2 <p, 



2 cos ip 



(33) 



(34) 



for the Co layers. 

In Fig.@]the layer-resolved effective anisotropy param- 
eters, K p s , are normalized to one Mn (layers no. 1-6) 
and Co (layers no. 7-12) atom, together with their one- 
site and two-site contributions. Reassuringly, when de- 
parting from the interface these functions show a stable 
convergence to the corresponding bulk values, K^ s =2.79 



and K P S 



K% 0i 

err 



K p £ s 

cn 



2.28 meV for IrMn 3 , 
for Co. In the vicinity 
of the interface, K P g S of Mn shows moderate oscillations 
and K P ff S of Co is still below the 0.01-0.02 meV range. 
The two-site anisotropy contribution at the IrMn3 layer 
nearest to the interface drops, however, to half of the 
bulk value. The absence of Mn-Mn interactions at the 
interface reduces the nearest neighbour coordination of 
the interface Mn spins from 8 to 6, the decrease observed 
indicates therefore that some softening of the intrinsic 
Mn anisotropy is found at the interface. Breaking of the 
cubic bulk symmetry, i.e., missing of the Co planes from 
one side, also results in an increase to about 0.22 meV of 
the Co K P g S at the interface. 



E. Dzyaloshinskii-Moriya interactions 

The antisymmetric components of the exchange in- 
teractions, Jfj = \ (Jij —Jjj), correspond to the well 
known Dzyaloshinskii-Moriya (DM) interaction! 31 ' 32 The 
DM interaction favours perpendicular alignment of mo- 
ments and for a given pair of spins, the antisymmetric 
component of the exchange interaction may be conve- 
niently represented by the DM vector, D^, 



e i^ij e j — Dij ■ (fii X Cj) 



(35) 



In the LI2 structure, the absence of inversion symme- 
try at the Mn lattice points gives rise to a DM inter- 
action with a magnitude of 1.24 meV between nearest 
neighbours, with Dij aligned along the (122) directions. 
Following the trend of the isotropic exchange interac- 
tions, the magnitudes of the DM interactions diminish 
very rapidly with distance; the third and fifth NN inter- 
actions in the bulk are 0.58 and 0.13 meV respectively. 
Clearly enough, in a non-collinear ground state of a bulk 
system, the asymmetric contribution to the exchange in- 
teraction might be non-zero only for different sublatticcs. 
However, due to symmetry, the inter-sublattice DM in- 
teractions also give rise to no overall magnetic anisotropy 
energy in bulk LI2 IrMn 3 . 

Calculations reveal that the symmetry breaking at the 
I1-M113/C0 interface gives rise to additional asymmetry 
in the exchange interactions at the interface. Although 
restricted to the atomic planes in the immediate vicin- 
ity of the interface, these interactions cause an anisotropy 
contribution for the reversal process discussed in the con- 
text of Fig. 2] that arises purely due to the presence of 
the interface and shows an angular dependence that is 
unidirectional, 



E™(<p)=E p (0) + K™coB<p 



(36) 



K^ M is plotted for the interface structure in Fig. [5] 
The reduced coordination of the interface Mn spins gives 
rise to a large anisotropic DM contribution to the system 
energy. Specifically, at the interface the antisymmetric 
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FIG. 5: (Colour online) Calculated layer-wise anisotropy con- 
stants due to DM interactions, Kp M , see Eq. (|36|1 . for the 
IrMn3/Co interface. 

exchange components give rise to an additional energy 
difference between the Tl and T2 Mn spin configurations, 
related to different chirality^ 

Concerning the Co layers, a value of —0.19 meV for 
the DM anisotropy energies (shown in Fig. [5]) denotes 
a unidirectional anisotropy favoring alignment of the Co 
spins perpendicular to, or away from, the interface plane. 
This anisotropy arises in the antisymmetric components 
of the interface Co-Mn interactions. The DM vectors for 
the interface exchange interactions are depicted as solid 
arrows in Fig. [6j where Dij is shown when i indexes a 
Co site and j a Mn site. 




FIG. 6: (Colour online) The arrows show the (111) in-plane 
projections of the DM vectors, Dij , for nearest neighbour in- 
terface Co-Mn interactions. The Tl configuration of the Mn 
spins is indicated by arrows within circles at the Mn posi- 
tions. Note that the plane of the Co atoms are above the 
IrMn3 layer. 

The interface DM vectors are found to lie almost en- 
tirely in the interface plane which favours a perpendic- 
ular spin alignment of the Co relative to the interface. 
Which of the two perpendicular directions is favoured 
depends upon the sense of the DM vectors and the Mn 



spin configuration: for the Mn configuration shown in 
Fig. [SJ Co alignment away from the IrMn3 layer is pre- 
ferred. Simultaneous inversion of the spins on all three 
Mn sublattices, which does not affect the chirality of the 
Mn spin structure, produces a configuration that is en- 
ergetically equivalent in the bulk, but which causes the 
preferred orientation of the Co layers at an interface to 
be reversed. Note, however, that such a change in the 
Mn spin configuration cannot be obtained by any single 
global rotation. 

IV. CONCLUSIONS 

We have introduced a new computational scheme to 
calculate magnetic anisotropy parameters and tensorial 
exchange interactions based on a Spin-Cluster Expansion 
and the Relativistic Disordered Local Moment picture 
of itinerant magnets. The method relies on the high- 
temperature paramagnetic phase, thus, it is, uniquely 
applicable to any system with arbitrary magnetic ground 
state. In particular, it is highly efficient in cases when the 
complex magnetic structure of the system is not known a 
priori from first principles total energy calculations, such 
as for heterogeneous interfaces and nanoparticlcs. More- 
over, the SCE-RDLM method can be used to calculate 
higher order multispin interactions. The main limitation 
of the method is that longitudinal spin-fluctuations are 
neglected: in case of large induced moments further ex- 
tension of the theory is needed^ 

We applied the new method to calculate magnetic 
anisotropics and spin-interactions at an IrMn3/Co(lll) 
interface. Our main observation is that the correspond- 
ing bulk properties are influenced in the immediate vicin- 
ity of the interface only. The exchange interaction be- 
tween the Mn and Co atoms are weak leaving the corre- 
sponding bulk magnetic structures fairly unaffected even 
at the interface. Symmetry breaking at the interface 
gives, however, rise to large anisotropy effects. Most 
prominently, the appearance of sizeable Dzyaloshinskii- 
Moriya interactions at the interface results in to a perpen- 
dicular coupling of the Co and Mn spins and to unidirec- 
tional exchange anisotropy that might substantially influ- 
ence the exchange bias in such systems^ - — By using the 
spin-model parameters derived in this work, we are going 
to present a detailed study to explore the spin-dynamics 
and the exchange bias properties of the IrMn3/Co(lll) 
interface. 
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Appendix A: Non-relativistic limit of the tensorial 
exchange interactions 

In order to derive a non-relativistic limit of the SCE- 
RDLM formalism for the two-site exchange interactions 
we first expand the logarithm to first order in Eq. ((26)) , 



the CPA condition (fT5l) can be rewritten as, 



LL 1 



with 



1 



Im 



2Li, L 



deTr (X iL t cij Xj,l< 



(Al) 



(A2) 



where, for convenience, but in this section only, we use 
complex spherical harmonics. In the non-relativistic case 
the transformation (|12p applies in spin space only, 



(A3) 



tt a &)=t i n6 ss ,+At i a a 



where and At, are now matrices in non-relativistic 
angular momentum (£, m)-space, a are the Pauli matrices 
and s =^, J, denotes the spin-index. Clearly, the usual 
spin-up and spin-down components of the i-matrices are, 



i = h - At, , 



(A4) 
(A5) 



(A12) 



The integration over angular variables can be explicitly 
performed in Eq. (|A2[) . 



(A13) 



X i>L = AX, a J d 2 e, Y* L (e* ) 



4-7r\ 3 

■y J AXj \25 L ,{i _i) a + - 2^,(1,1) <t- + £l,(i,o) &z) > 

with cr± = <t x ± iijy. The trace in Eq. (| Al|) can then 
be recast into a trace in angular momentum space, Tr^, 
and a trace in spin space, Tr s , 



Tr (X_,i. L,.,,A\,j/L, .,,) 

— ' ' ' (AX ilciij AA ;Z , „ ) 

Tr s [(2^,(1,-1) ff+ - 2( 5l,(i,i) o-- + <*l,(i,o) °o) 

( 2 ^L',(i,-i) o"+ - 2(5 L / : ( ia ) cr_ + (5z,',(i,o) fo)] ■ 

(A14) 



while for the paramagnetic case the effective i-matrices 
and SPO matrices become diagonal in spin space, 



tl 



—ij,c u ss 



(A6) 
(A7) 



Since also the excess scattering matrices take the form of 
JH and EH) 



By using the algebra of the Pauli matrices, the trace 
in spin space can easily be performed and the two-site 
terms in the spin-Hamiltonian will be transformed into a 
pure isotropic Heisenberg interaction, 



LL' 



(A15) 



AT (el) 



X_ i 5 SS > 



AX, ~3 , 



(A8) 



with 



with 



Ki — 2Lm + AXi , 
Xj = X i0 -AX i , 



and 



X s , = 



fa-®)- 1 ) 1 -!* 



-i -l 



(A9) 
(A10) 



(s=t4) 



(All) 



Jii = 



- 1 Im de Tr L (Al, r, ^ AX^ r^) (A16) 



- — Im / de TrL 

Z7T 



A,- — X_i j Z c ,i 3 



(A17) 



The above formula is widely used in non-relativistic cal- 
culations of the exchange interactions. 
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